Interplay of electron-phonon interaction and strong correlations: DMFT+E approach 
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We discuss interaction of strongly correlated electrons (described within the Hubbard model solved by dynamical mean-field 
theory (DMFT)) with Debye and Einstein phonons using recently developed DMFT+X computational scheme. Electron-phonon 
interaction (EPI) is analyzed in adiabatic approximation (assuming the validity of Migdal theorem), allowing the neglect of EPI 
vertex corrections. This approach is valid for EPI coupling constant A < Sf/copi, ~ 10, where Sf is Fermi energy and (±> p h is Debye 
or Einstein frequency. For moderate values of A only small changes in the electronic density of states are observed in DMFT+S 
approximation for both weakly and strongly correlated metallic regimes. Metal-insulator (Mott) transition due to the increase of 
Hubbard interaction U is slightly inhibited by EPI. Our main aim is to discuss the interplay of "kinks" in electronic dispersion due to 
EPI and recently discovered kinks of electronic origin. For the certain region of model parameters coexistence of phonon "kinks" in 
electronic dispersion with purely electronic "kinks" is readily observed and we formulate some simple criteria of such coexistence. 
However, for most general combinations of model parameters phonon "kinks" make electronic "kinks" hardly observable. In the 
general case an increase of Hubbard interaction U rapidly suppresses the slope of electronic dispersion within the phonon "kink". 
These results are important for deeper understanding of the shape and evolution of electronic dispersions in strongly correlated 
systems such as copper oxides, where different kinds of "kinks" were recently observed in ARPES experiments. 
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1. Introduction 

The problem of the interplay of strong electronic correlations 
with electron-phonon interaction is of central importance in the 
physics of highly correlated systems. Actually there is rather 
long history of such studies, e.g. one of the most popular mod- 
els for electron-phonon interaction (EPI) in strongly correlated 
systems is the so-called Hubbard-Holstein model (HHM). The 
Hubbard model [1] itself describes local Coulomb interaction 
of electrons on a lattice including e.g. Mott-Hubbard metal- 
insulator transition. On the other hand Holstein model contains 
local linear displacement-to-density interaction of conducting 
electrons with local (Einstein) phonon modes yfl. 

Active investigations of the properties of the HHM were 
undertaken in the framework of dynamical mean-field theory 
(DMFT) [3], which is non-perturbative approach with respect 
to interaction parameters of the Hubbard model. Among many 
others one should mention DMFT solution of HHM for the 
case where impurity solver used was the numerical renormal- 
ization group (NRG) |4]. The mapping of HHM to Anderson- 
Holstein impurity was first performed by Hewson and Mayer 
iH. It was shown that using NRG one can compute in a numer- 
ically exact manner total electron-phonon contribution to the 
self-energy of the problem, thus making solution of the HHM 



non-perturbative also with respect to electron-phonon coupling 
strength. 

However, up to now there are apparently no studies of 
strongly correlated electrons interacting with Debye phonons. 
It is even more surprising in view of the widely discussed 
physics of kinks in electronic dispersion observed in ARPES 
experiments 40-70 meV below the Fermi level of high- 
temperature superconductors |6], which are often attributed to 
EPI 1 7]. To our knowledge problem of kink formation on elec- 
tronic dispersion caused by EPI in strongly correlated systems 
was briefly discussed within HHM in papers by Hague [9] and 
Koller etal. fl. 

In this paper we consider the influence of Debye or Einstein 
phonons on the weakly and strongly correlated electrons within 
our recently developed DMFT+2 approach, studying electron 
dispersion and density of states (DOS), in particular close to 
Mott-Hubbard metal insulator transition. We analyze in de- 
tails how EPI affects electronic dispersions in correlated metal 
and discuss the interplay of recently discovered kinks of purely 
electronic nature in electronic dispersion IU0I1 and usual phonon 
kinks in the electronic spectra. 



2. DMFT+E computational details 
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The major assumption of our DMFT+Z approach is that the 
lattice and time Fourier transform of the single-particle Green 
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function can be written as: 



Eq. (O can be rewritten as (cf. similar analysis in Ref. 118]) 
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(1) 



where e(p) is the bare electron dispersion, is the lo- 
cal self-energy of DMFT, while E p (s) is some "external" (in 
general case momentum dependent) self-energy. Advantage 
of our generalized approach is the additive form of the self- 
energy (neglect of interference) in Eq. ([T) ifTll IT2L [l3ll . It al- 
lows one to keep the set of self-consistent equations of stan- 
dart DMFT |3j]. However there are two distinctions. First, 
on each DMFT iterations we recalculate corresponding "ex- 
ternal" self-energy "Z p (fi,s, [2(e)]) within some (approximate) 
scheme, taking into account interactions e.g. with collective 
modes (phonons, magnons etc.) or some order parameter fluc- 
tuations. Second, the local Green's function of effective impu- 
rity problem is defined as 
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N ^ e + A i-e(p)-Z(e)-Sp(e)' 
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at each step of the standard DMFT procedure. 

Eventually, we get the desired Green function in the form of 
(Q3, where £(e) and S p (e) are those appearing at the end of our 
iteration procedure. 

To treat electron-phonon interaction for strongly correlated 
system we just introduce 2 p (e) = 2 p ;,(e, p) due to electron- 
phonon interaction within the usual Frohlich model. To solve 
single impurity Anderson problem we use NRGjj]. All calcu- 
lations are done at nearly zero temperature and at half filling. 
For "bare" electrons we assume semielliptic DOS with half- 
bandwidth D. 

According to the Migdal theorem in adiabatic approxima- 
tion 11411 we can restrict ourselves with the simplest first or- 
der contribution to E p /,(e, p), neglecting vertex corrections due 
electron-phonon coupling which are small over adiabatic pa- 
rameter (a>D, coq)/sf ^ 1 [14]: 
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e + id + li - e(p + k) - £(s + co) - S p /,(e + co, p + k) 
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where g is the usual electron-phonon interaction constant, 
wo(k) is phonon dispersion, which in our case is taken as in 
the standard Debye or Einstein model 
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Here u is the sound velocity, cod, coq are Debye and Einstein 
frequencies with cut-off ko of the order of Fermi momentum 
Pf- 

Actually 2 ; ,;,(e, p) defined by Eq. © has weak momentum 
dependence which we can omit and continue only with signif- 
icant frequency dependence. For Debye spectrum of phonons 
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with a characteristic frequency co c - p F u of the order of cod, 
while for Einstein spectrum: 
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where E s = e - E(e) - £ p /,(e). For the case of semielliptic 
non-interacting DOS No(s) with half -bandwidth D we get: 



7(e) = ^{E s - J E* - D 2 ), 



(8) 



It is convenient to introduce the dimensionless electron-phonon 
coupling constant as lfl8ll : 



a) 2 k 2 
A D = g 2 N Q (s F )-^, A E = g 2 N Q (e F )-\ 



(9) 



To simplify our analysis we shall not perform fully self- 
consistent calculations neglecting phonon renormalization due 
to EPI B18I1 . assuming that the phonon spectrum <(4j is fixed by 
the experiment. 

3. Results and discussion 

Let us start from comparison between pure DMFT 
and DMFT+2,,/, DOSes for strong (U/2D=1.25) and weak 
(U/2D=0.625) Hubbard interaction presented in Fig.[T]on upper 
and low panels correspondingly. Dimensionless EPI constant 
(O used in these calculations was Ad = /l£=0.8, while Debye 
and Einstein frequencies were taken to be wd=«o=0.125D. In 
both cases we observe some spectral weight redistribution due 
to EPI. For U/2D=1.25 (upper panel of Fig. [TJ we see the well 
developed three peak structure typical for strongly correlated 
metals. In the energy interval +cod, too around the Fermi energy 
(which is taken as zero energy at all figures below) there is al- 
most no difference in the DOS quasiparticle peak line shape ob- 
tained from pure DMFT and DMFT+Z,,/,. However outside this 
interval DMFT+2^/, quasiparticle peak becomes significantly 
broader with spectral weight coming from Hubbard bands and 
it is more pronounced for the case of Einstein phonons. This 
broadening of DMFT+I^/, quasiparticle peak leads as we show 
below to inhibiting of metal to insulator transition. In the case 
of U/2D=0.625 there are no clear Hubbard bands formed but 
only some "side wings" are observed. Spectral weight redis- 
tribution on the lower panel of Fig. Q]is not dramatic, though 
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Energy in units of D 
Figure 1: Comparison of DOSes obtained within standard DMFT (dashed 
lines) and DMFT+Sp/, methods (solid lines) for Debye (black) and Einstein 
(red) phonons for strong (upper panel, U/2D=1.25) and weak (lower panel, 
U/2D =0.625) Hubbard interaction regimes. Dimensionless electron-phonon 
coupling constant Ap = Ae=Q.&. 



qualitatively different from the case of U/2D=1.25. Namely, 
main deviations between pure DMFT and DMFT+Sp/, happen 
in the interval ±o>d, where one can observe kind of "cap" in 
DMFT+Zp/, DOS. Corresponding spectral weight goes to the 
energies around ±U, where Hubbard bands are supposed to 
form. The lineshape of the "cap" is slightly different for the 
Debye and Einstein phonons due to different behavior of S p /, 
(Eqs. (0, (0) at energies ±u>d, ojq- For Einstein phonons lm£ p i, 
at these energies sharply drops down to zero. This leads to sharp 
cusps of DOS at ±ojq as shown at the insert in Fig.Q] 

In Fig. |2] we compare the behavior of pure DMFT and 
DMFT+Sp/, DOSes for different U/2D values close to Mott- 
Hubbard metal-insulator transition for the case of Debye 
phonon spectrum. For U/2D=1.56 both standard DMFT and 
DMFT+Z^z, produce insulating solution. However there is 
some difference between these solutions. The DMFT+S p a Hub- 
bard bands are lower and broader than DMFT ones because 
of additional interaction (EPI) included. With decrease of U 
for U/2D=1.51 and 1.47 we observe that DMFT+2,,/, results 
correspond to metallic state (with narrow quasiparticle peak at 
the Fermi level), while conventional DMFT still produces in- 
sulating solution. Only around U/2D=1.43 both DMFT and 
DMFT+Z^, results turn out to be metallic. Overall DOSes line- 
shape is the same as discussed above. These results show that 
with the increase of U finite EPI slightly inhibits Mott-Hubbard 
transition from metallic to insulating phase. For the case of 
Einstein phonons the MIT is inhibited even stronger. This re- 
sult is similar to what was observed for the HHM in weak EPI 
regmeJnlEII!. For more deep insight into these results on 
DOS we have also analyzed the fine structure of corresponding 
self-energies 2(e) and Epft(s). Relevant details can be found in 
Ref. dl. 

Now we address the issue of a sudden change of the slope of 
electronic dispersion, the so-called kinks. It is well known that 
interaction of electrons with some bosonic mode always pro- 
duces such a kink. In the case of EPI typical kink energy is just 
the Debye or Einstein ojq frequency. Kinks of purely elec- 
tronic nature were recently reported in Ref. jlOll. The energy 
of purely electronic kink as derived in Ref. [ 10] for semielliptic 
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Figure 2: Sequence of DOSes obtained within standard DMFT (dashed lines) 
and DMFT+X,,/, for Debye model (solid lines) methods close to metal-insulator 
transition (from top-left to bottom right) with Ap=Q.8. 



bare DOS is given by 

u* =Z FL (V2- 1)D, 



(10) 



where D is the half of the bare bandwidth of electrons and 
Zfl = (1 - | s=e ) ' lS Fermi liquid quasiparticle weight. 
The rough estimate of of is given by the half-width of quasi- 
particle peak of DOS at its half-height. 

Our calculations clearly demonstrate that electronic kinks are 
hardly observable on the background of phonon kinks and spe- 
cial care should be taken to separate them by rather fine tuning 
of the parameters of our model. To clarify this situation we 
introduce an additional characteristic of the kink — the shift 
of electron dispersion in momentum space dp at kink energy. 
From simple geometry we estimate for phonon kinks 
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where v F is the bare Fermi velocity and Ad,e was defined in 
Eq. (0. For electronic kink the similar estimate is 
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where Zo is quasiparticle weight in the case of absence of elec- 
tronic kinks (the same as Z cp defined in Ref. lfl(ilo . Velocity v* F 
is the Fermi velocity of initial dispersion, but it can not be just 
a bare one. As was reported in Ref. II 1(111 electronic kinks can be 
observed only for rather strong Hubbard interaction when three 
peak structure in the DOS is well developed and electronic dis- 
persion is strongly renormalized by correlation effects. This 
renormalization is determined by A e defined in Eq. (fTTT i. which 
can be seen as kind of dimensionless interaction constant. In 
the case when both slopes on the Fermi level and out of ±cf 
energy interval are equal there will be no electronic kink at all. 

Now we can choose parameters of our model to make both 
kinks simultaneously visible. First of all one should take care 
that co D <k co*. For U/2D=1 with U=3.5 eV we get of ~ 
0.1D and a reasonable value of Debye (or Einstein) frequency 
is a»o(or ojq) ~ 0.01D. To make phonon kink pronounced at 
such relatively low Debye (or Einstein) frequency (cf. Eq. (fTTb ) 
we have to increase EPI constant and we take Ad = Ae=2.0. To 
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Figure 3: (Color online) Quasiparticle dispersions obtained from standard 
DMFT (black dashed line) and DMFT+S pA for Debye (solid black line) and 
Einstein (solid red line) phonons along the part of high symmetry direction 
T - (n,7T,n). 

demonstrate coexistence of both these types of kinks we plot the 
energy dispersion of simple cubic lattice with nearest neighbors 
transfers only, along the high symmetry direction V - {n, n, n) 
lioll . In Fig. [3] we show dispersion along this direction close to 
the Fermi level. The difference of lineshapes of Debye and Ein- 
stein kinks is illustrated at the insert in the Fig. [3] As discussed 
above "Einstein" kink is more sharp. 

Finally we address to the behavior of phonon kinks in elec- 
tronic spectrum as function of Hubbard interaction U. As U/2D 
ratio grows Fermi velocity in Eq. (fTTb goes down, so that mo- 
mentum shift of kink position 8p moves away from pf , while 
kink energy remains at loq- This is confirmed by our direct 
DMFT+S^/, calculations producing the overall picture of spec- 
trum evolution shown in Fig. [4] 
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Figure 4: (Color online) Quasiparticle dispersions around Fermi level with De- 
bye phonon kinks obtained from DMFT+Z,,/, calculations for different interac- 
tion strengths U/2D= 0.5, 0.75, 1.0; A D = 0.8, w D =0.1D. 



4. Conclusion 

This work is a first attempt to analyze strongly correlated 
electrons, treated within DMFT approach to the Hubbard 
model, interacting with either Debye or Einstein phonons. EPI 
was treated within the simplest (Migdal theorem) approach in 
adiabatic approximation, allowing the neglect of vertex correc- 
tions. DMFT+Sp/j approach allows us to use the standard mo- 



mentum space representation for phonon self-energy (O, while 
the general structure of DMFT equations remains intact. 

Mild EPI leads to rather insignificant changes of electron 
density of states, both in correlated metal and in Mott-insulator 
state, slightly inhibiting metal to insulator transition with in- 
crease of U. However, kinks in the electronic dispersion due 
to EPI dominate for the most typical values of the model pa- 
rameters, making kinks of purely electronic nature, predicted 
in Ref. [10], hardly observable. Special care (fine tuning) of 
model parameters is needed to separate these anomalies in elec- 
tronic dispersion in strongly correlated systems. We have also 
studied phonon kinks evolution with the strength of electronic 
correlations demonstrating the significant drop in the slope of 
electronic dispersion close to the Fermi level with the growth 
of Hubbard interaction U. Quantitative difference of results for 
the cases of Debye and Einstein phonon spectra was observed 
both in DOS and kink behavior in electronic dispersion. 
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